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Electro-optical properties of deformed helix ferroelectric liquid crystal (DHFLC) 
cells are studied by using a general theoretical approach to polarization gratings in 
which the transmission and reflection matrices of diffraction orders are explicitly 
related to the evolution operator of equations for the Floquet harmonics. In the 
short-pitch approximation, a DHFLC cell is shown to be optically equivalent to 
a uniformly anisotropic biaxial layer where one of the optical axes is normal to 
the bounding surfaces. For in-plane anisotropy, orientation of the optical axes and 
birefringence are both determined by the voltage applied across the cell and represent 
the parameters that govern the transmittance of normally incident light passing 
through crossed polarizers. We calculate the transmittance as a function of the 
applied voltage and compare the computed curves with the experimental data. The 
theoretical and experimental results are found to be in good agreement. 
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I. INTRODUCTION 



A polarization grating (PG) can generally be described as an optically anisotropic layer 
characterized by the anisotropy parameters that periodically vary in space along a line in the 
plane of its input face. Unlike conventional phase and amplitude diffraction gratings, PGs 
act by locally modifying the polarization state of light waves passing through them. Ow- 
ing to the one-dimensional (ID) in-plane periodicity, this introduces periodically modulated 
changes of the polarization characteristics giving rise to polarization-dependent diffraction. 
In particular, the latter implies that a PG divides a monochromatic plane wave into differ- 
ently polarized diffracted waves. 
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Over the past decade PGs have been attracted much attention due to a unique combi- 
nation of their optical properties: (a) it is possible to achieve 100% diffraction into a single 
order; (b) diffraction efficiencies are highly sensitive to the incident light polarization; and 
(c) the state of polarization of diffracted orders is determined solely by the parameters of 
a PG 0-0] ■ There are numerous applications in a variety of fields, including polarimeters, 
displays, polarizing beam splitters, beam steering and polarization multiplexers where PGs 
have been found to be useful (for a recent review, see the article jEf and the monograph 0). 

There are different technologies to fabricate PGs. For example, computer-generated 
subwavelength-period metal-stripe gratings with spatially periodic fringe orientation [tJ and 
space- variant dielectric subwavelength gratings formed by discrete orientation of local sub- 
wavelength grooves jl, @] are produced using advanced photolithographic and etchingtech- 



niques. Such gratings were employed to perform real-time polarization measurements |8MlO 
They are also used to demonstrate polarization Talbot self-imaging and Pancharatnam- 
Berry phase optical elements (l~2l . fl3j | . 

Polarization holography provides another well-known method to produce PGs |6j. It uses 
two differently polarized light beams to record the spatially modulated polarization state of 
the resultant light field on suitable media such as azobenzene containing polymer systems 
and silver-halide materials. 

The holographic technique has been extensively used to create polarization gratings in 
liquid crystal (LC) cells with photosensitive aligning substrates such as linear photopoly 



merizable polymer layers [H], 15], azo-dye films (l6l |. azo-dye doped polyimide , and 



azobenzene side-chain polymer layers [19]. 

In this method, irradiation of the substrate with a holographically generated polariza- 
tion interference pattern gives rise to spatially modulated light induced ordering in the 
photoaligning layer. This ordering manifests itself in the effect of photoinduced optical 
anisotropy and determines the anchoring properties of the layer such as its (polar and 
azimuthal) anchoring strengths and the easy axis orientation (see, e.g., Refs BS and 
references therein). 

The anchoring parameters of the photoaligning film thus undergo periodic variations 
across the substrate face leading to the formation of orientational structures in the liquid 



crystal cells 23] characterized by spatially periodic distributions of the liquid crystal director, 
d = (d x ,d y ,d z ), which is a unit vector that defines a local direction of the preferential 
orientation of LC molecules. In liquid crystals, the elements of the dielectric tensor, e, can 



be expressed in terms of the LC director [2j 

= e ± ((% + u a didj), u a = (ey - ej_)/e ± , (1) 

where 8ij is the Kronecker symbol, u a is the anisotropy parameter and n_|_ = n Q = ^fie± 
(ny = n e = ^JJie~\\) is ordinary (extraordinary) refractive index (the magnetic tensor of LC 
is assumed to be isotropic with the magnetic permittivity n). So, the periodic orientational 
LC configurations define the so-called liquid crystal polarization gratings (LCPG). These 
gratings will be of our primary interest. 

Common methods most generally employed to derive theoretical results for PGs typically 
rely on on the well-known Jones matrix formalism and its modifications 1^]. These 



results are limited by their assumptions to large gratings periods and normal incidence. In 
addition, using Jones calculus implies neglecting multiple reflections. 

In this work we present the theoretical approach to PGs that can be regarded as a 



generalized version of our method developed in Refs. [25l . |26| for stratified anisotropic media 



3 



and goes beyond the limitations of Jones calculus. We apply the method for systematic 
treatment of the technologically important case of the deformed helix ferroelectric (DHF) 



liquid crystals |27|, |28| where the director rotates about a uniform twist axis parallel to the 



substrates forming the ferroelectric LC (FLC) helical structure (see Fig. [I]). 



* 
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Figure 1. Helical structure of deformed helix ferroelectric liquid crystals. The director rotates 
around the helix axis (the x axis) with cone angle 9t- 

In DHFLC cells, the FLC helix is characterized by a short submicron helix pitch, P < 
1 /mi, and a relatively large tilt angle, 9 t > 30°. Note that, in the case of surface stabilized 
FLC cells, the helix pitch of a FLC mixture is typically greater than the cell thickness, so 
that the bulk chiral helix turned out to be unwound (suppressed) by the boundary conditions 
at the substrates j29[. By contrast to this, a DHFLC helix pitch is 5-10 times smaller than 
the thickness. This allows the helix to be retained within the cell boundaries. 

Electro-optical response of DHFLC cells exhibits a number of peculiarities that make 
them useful for LC devices such as high speed spatial light modulators [3(J [3l[ and colour- 
sequential liquid crystal displays [32]. So, in this study, our goal is to examine electro- 
optical properties of DHF liquid crystals based on the general theoretical approach describing 
polarization gratings. 

The layout of the paper is as follows. 

In Sec. [Ill we begin with Maxwell's equations for the lateral components of the electric 
and magnetic fields and derive a set of equations for the Floquet harmonics representing 
diffracted waves. The relations linking the transmission and reflection matrices of diffraction 
orders and the evolution operator of the system for the harmonics are deduced in Sec. Ill CI 

Electro-optical properties of DHFLC cells with the subwavelength helix pitch are studied 
in Sec. IIHI Experimental details are given in Sec. IIII Al where we describe the samples and 
the setup employed to perform measurements. In Sec. IIII 13^ the general theory of Sec. HTlis 
used to examine how the optical anisotropy parameters and the transmission coefficients of 
the DHFLC cells depend on the applied electric field. In particular, it is found that in the 
short-pitch approximation the DHFLC polarization gratings can be represented by uniformly 
anisotropic biaxial layers. For the electric field dependence of the light transmittance through 
the cell placed between crossed polarizers, the results of electro-optic measurements are 
compared with the theoretically computed curves in Sec. IIII CI 

Finally, in Sec. \1V\ we present the results and make some concluding remarks. Technical 
details on derivation of Maxwell's equations for the lateral (in-plane) components of elec- 
tromagnetic field are relegated to Appendix Theory of stratified media is recapitulated 
in Appendix IBl 
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II. THEORY OF POLARIZATION GRATINGS 



In this section we shall generalize the theoretical approach developed in Refs. [25|, [26| so 
as to treat the light transmission problem for a polarization grating in the slab geometry 
illustrated in Fig. [2j In this geometry, as is indicated in Fig. [21 the z axis is normal to the 
bounding surfaces of the layer: z = and z = D, the grating with the grating pitch, A g , and 
the grating wave vector, 

kg = kgX, kg = —, (2) 

where k g is the grating wave number, is characterized by the condition of in-plane periodicity 
for the elements of the dielectric tensor, e: 



[x + A g ) = e ij (x). (3) 



Expression for the grating wave vector §2§ defines the x-z plane as the plane of grating. 

Throughout this paper we deal with harmonic electromagnetic fields characterized by 
the frequency, u, (time-dependent factor is exp{— iut}) and the free-space wave number, 
k vac = oj/c. So, the starting point of our theoretical considerations is the Maxwell equations 
for a harmonic electromagnetic wave written in the form: 

V x E = ifik VSiC H, (4a) 

V x H = -^ vac D, (4b) 



where D = e ■ E is the electric displacement field. 
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We shall also assume (see Fig. [2]) that the medium surrounding the layer is optically 
isotropic with the dielectric constant e m , the magnetic permittivity /i m and the refractive 
index n m = ^fi m e m . So, for a plane wave travelings along the wave vector 

k = k m k = k z z + k p , k p = k p [cos(0) x + sin(0) y] , k m = n m k vac , (5) 
the electromagnetic field {E, H} is given by 

{E, H} = {E(k), H(k)} exp[*(k ■ r)], (6a) 
E(k) = E p e x (k)+E s e y (k), (6b) 

^ H(k) = k x E(k) = E p e y (k) - E s e a (k), (6c) 
n m 

where the unit vectors k = (sin 9 cos 0, sin 9 sin 0, cos 9) , e x (k) = (cos 9 cos 0, cos 9 sin 0, — sin 9) 
and G y {k) = (— sin0, cos0, 0) expressed in terms of the polar (9) and azimuthal (0) angles 
form an orthogonal basis. 

The incoming incident wave {E inc , H inc } is represented by a plane wave (Q, {E, H} = 
{E inc , H inc }, propagating along the wave vector 

k = kj nc = k\, -*z + k p , kS z ^ = yj k^ — k^ (7) 

in the half space z < bounded by the input face of the grating. In this case, the polar 
angle, < 9 = 9- mc < 7r/2, is the angle of incidence, whereas the azimuthal angle, < = 
inc < 27r, is the angle between the grating plane and the plane of incidence. 



A. Maxwell's equations for lateral components 



Now we write down the representation for the electric and magnetic fields, E and H, 

E = E z z + E P , H = H z z + z x H P , (8) 

where the components directed along the normal to the bounding surface (the z axis) are 
separated from the tangential (lateral) ones. In this representation, the vectors Ep = 

E x yL + Eyf = yjijj anc ^ Hp = Hxz = ^ ^ J are parallel to the substrates and give the 

lateral components of the electromagnetic field. Similar decomposition for the differential 
operator that enter the Maxwell equations (j3J) is given by 

Km V = z d T + i V p , V; = zx V p> (9) 

where r = k vac z; V p = -i k^ c (±d x + yd y ) = (V x , V v ) and = (V£, V^) = (-V„, V x ). 

We can now substitute Eqs. jS]) and (jSJ) into the system (jlj) and follow the algebraic pro- 
cedure described in Appendix |X] to eliminate the z-components of the electric and magnetic 
fields. As a result, the ^-components, E z and H z , turned out to be expressed in terms of the 
lateral components (see Eqs. ( ]A5j) and (IA6p ) . The resultant system of differential equations 
for the lateral components, Ep and Hp, can be written in the following matrix form: 
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where the elements of the matrix differential operators Aiij are given by 

M% ] = -V a -[e;le z pi M^f = ^ - V Q • e" 1 • V* (lfa) 

= -ea.£Vp, Mtf^e^-^Vi-V^ (ffb) 

and the elements of the effective dielectric tensor f]A8|) that enter the operator are 

e< 3 = e «/3 - e az^e zf3 , a,Pe {x, y}. (12) 



B. Floquet harmonics 

From Eq. ([3]), the dielectric tensor is a periodic function of x, so that it is represented by 
its Fourier series expansion 

oo oo 

£ = £ nexp[in(k g ■ r p )) = ^ e n exp(ink^c), (13) 

n=— oo n=— oo 

where r p = (x,y,0). Similar remark applies to the coefficients that enter the differential 
operators (fTT|) . 

The representation of Floquet harmonics [33[ for solutions of the system (flQj) 



F(r) = F(r p , r) = ^ F «( r ) exp{i(k n • r p )}, k n = k p + nk g = fc vac q n , (14) 

n=—oo 

q n = g n (cos0 n x + sin0 n y) = (q^\ q y n \ 0), (15) 

where k p is the tangential component of the wave vector of the incident plane wave defined 
in Eq. ([7]), is another consequence of the periodicity condition ([3]). 

On substituting the expansion over the Floquet harmonics ( fT4"|) into the system ( TTOj) . we 
derive a set of matrix equations for the Floquet harmonics 

oo 

-id T ¥ n (r)= M nm (r)-F m (r), (16) 



m=— oo 



/M (11) M (12) \ 1 f As — 
Mnm = tv/t?21) A/r ?ii) = T~ / ex P{-^( k " • r )> [-^ exp{i(k m • r)}]dx, (17) 

\ lVl nm lVl nm / iV g Jo 

where the 2x2 block matrices Mnm are given by 

[m£J] a , = -^ } /?S- m) ' I M 2] = -/e m) 4 m) , ( 18a ) 

^A m -^ n) ^r m) ^ m) ' (18b) 
e iV m) ~ ^^nm^U™', p„ = Z X q n , (18c) 



I a j3 ~ e ce/3 ^ °nm Va P/3 

and rjiz , (3^ and e^l are the Fourier coefficients for e~£, t~}e^ and ef), respectively 
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General solution of the system ([175]) 

oo 

F -( r )= U nm (r,r )-F m (r ) (19) 



m=— oo 



can be conveniently expressed in terms of the evolution operator defined as the matrix 
solution of the initial value problem 

oo 

-id T \J nm (r,T ) = M nfc (r) ■ U fem (r, r ), (20a) 

k=— oo 

U nm (r , r ) = I 4 S nm , (20b) 

where I n is the n x n identity matrix. 

The Floquet harmonics, F n , represent the electromagnetic field of the diffracted waves 
with the integer nfZ giving the diffraction order. From the system f JT6|) . it can be inferred 
that the effect of inter-harmonics coupling responsible for diffraction is solely caused by the 
periodic modulation of the dielectric tensor ffT3]) . 

In the ambient medium with = e m 5ij and \i = /x m , the harmonics are decoupled and, 
for the nth diffraction order, represent plane waves propagating along the wave vectors with 
the tangential component (IT5|) . The Floquet harmonics of such waves, is given by 

F M frl _ v (q) ( exp{?Q m (g n ) r} \( 

n [ ' ~ m[qn) V exp{-*Q m (g n )r} ) E W ' 1 ' 



Qm(gn) = qm(q n )h, = V n li ~ 9n> ( 22 ) 

where V m (q n ) is the eigenvector matrix for the ambient medium given by 



RtUr,) WE m -cr 3 E m 



V m (q„) = T„ t («V m(9 „) = {"J" Rt J M j ^« — j , (23) 

E ™-{ lj' ^^-^0 g m (g„)J' (24) 

Rt(0)=f COS ^- sin /V (25) 
v ' \ sm0 cos0 y 

and {<Ti, <r 2 , cr 3 } are the Pauli matrices 

i\ A") -/A (\ o 

1 o 



»> = U o ■"■= o-i ■ (26 » 



From Eq. ( I2ip . the vector amplitudes and E^ 1 correspond to the forward and back- 
ward eigenwaves with ki + ^ = +k vs , c q m and ki = —k vac q m , respectively. In the half space 
z < 0, these describe the incident and reflected waves 



^i(inc) 

E+ | z < = E;™ c = 5„ E inc = 5 n0 [ J? inc) , (27) 



EL%<o = E2=(^ fl) ). (28) 

-L-/S, n 
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Clearly, equation ( 12 7p implies that the incident wave is the forward eigenwave of the zeroth 
order. 

In the half space z > D after the exit face of the grating, the only wave of the nth order 
is the transmitted plane wave 

t p (trm)\ 
^(trm) I > \z>D - U. {/V) 

Note, that, at sufficiently large diffraction order with q n > n m and q m = i\q m \ (see Eq. (|22|) ). 
the reflected and transmitted waves become evanescent. In this case, the ^-components of 
the wave vectors, = k v&c (-q m z + q n ) and = C(9 m z + q„), are imaginary. 



C. Computational procedure 

We can now define the transmission and reflection matrices through the linear input- 
output relations 

-^trm = ' E; nc , ^refl = ^ n ' ^inc (30) 

linking the nthe order for the transmitted and reflected waves to the incident wave and, 
following the line of reasoning presented in Refs. 25|, |26), relate these matrices and the evo- 
lution operator given by Eq. (|20|) . To this end, we use the boundary conditions requiring the 
tangential components of the electric and magnetic fields to be continuous at the boundary 
surfaces: F n (0) = Fl m) (0 - 0) and F n (h) = F^ih + 0), and apply the relation flU]) to the 
anisotropic layer of the thickness D to yield the following result 

F^\h + 0) = Unk(h, 0) • F<: m) (0 - 0), h = k vac D. (31) 

k 

By using Eq. ( 12 ip . equation (13 ip can be conveniently recast into the form 

E- (n) \ fE (k) 



T7i(n) 



E w «*( o j' (32) 



where W n ^ is given by 



W (H) W (12) 



Wnk = [V m (q n )]^ • [\J-\h, 0)U • V m (q fe ) = ? 2 % # 2) . (33) 

\ nk nk J 

and defines the linking matrix. 

At n 0, the relation (13"2"|) gives the system for the transmitted wave orders 

-W^.E£ = ^W^.Eil n^0 (34) 

k^0 

By solving this system we obtain the diffracted waves 

Etrm — T n • E|°^ (35) 
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linearly related to the zeroth order (non-diffracted wave) through the efficiency matrices T n . 
Obviously, Eq. ([35]) implies that, by definition, T = 12- 
From Eq. f )32|) at n = 0, we obtain the relation 



(ii) 

Ok 



E 



(36) 



linking the vector amplitudes of the incident and transmitted waves. Substitution of the 
matrices from Eq. f l33|) into Eq. f[3"6"l) provides the transmission matrix (I3"U1) in the following 
form: 

-i 



(ii) 

Ofc 



(37) 



From the relation (|32|) applied to the case of reflected waves, we derive the reflection matrix 



(38) 



expressed in terms of the transmission matrices ([2 

In the limiting case of uncoupled harmonics, where W n 



recover the results for stratified media obtained in Refs. [25|, [26j. It will suffice to note that, 



^nfeW nn , it is not difficult to 



at W 



(ii) 

On 



0, the system 



gives the efficiency matrices T n = <5 n0 12 that point to the 



absence of diffracted waves. 




Figure 3. Geometry of a deformed helix FLC cell. 



III. DEFORMED HELIX FERROELECTRIC LIQUID CRYSTAL CELLS 



In this section we present the experimental results on the transmittance of light passing 
through crossed polarizers measured as a function of the applied electric field (voltage) in 
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DHFLC cells. In order to interpret the experimental data, the theoretical approach of Sec. [TT] 
is applied to chiral smectic helical configurations with the subwavelength helix pitch. 



A. Experiment 



1. Sample preparation 



The photo-alignment technique described in Ref. [34| was used for producing the FLC 
cells. Geometry of the cells is schematically depicted in Fig [3J We used the cells with the 
size of 13 x 13 mm 2 , the thickness of the glass substrate 1.1 mm, electrodes area 10 x 10 mm 2 , 
and the cell thickness (gap) 50 //m and 130 /xm. 

Following the method of Ref. 34J , ITO surfaces of FLC cells were covered with a 10-20 nm 
photo-aligning substance - azobenzene sulfonic dye SD-1 layers. The azo-dye solution was 
spin-coated onto ITO electrode and dried at 155 °C. 

The surface of the coated film was illuminated with linearly polarized UV light using a 
super-high-pressure Hg lamp through an interference filter at the wavelength 365 nm and 
a polarizing filter. The intensity of light irradiated normally on the film surface during 30 
minutes was 6 mW/cm 2 . 

In our experiments we used the FLC mixture FLC-576A (from P. N. Lebedev Physical 
Institute of Russian Academy of Sciences) as a material for the DHFLC layer. The mixture 
was injected into the cells in the isotropic phase by capillary action. This mixture has the 
helix pitch P m 300 — 330 nm. Hence, in visible spectral range, the light scattering turns 
out to be completely suppressed [35[ except that the applied voltage is close to the critical 
voltage of the helix unwinding. So, such mixtures exhibit pronounced effects of electrically 
controlled phase retardation and birefringence An eS (V) that can be clearly detected. 




Figure 4. Experimental set-up for electro-optical measurements of DHFLC cells. 



2. Experimental setup 

The measurements of the light transmittance in relation to the applied voltage were 
performed in the automatic regime. For this purpose, the measurement complex device was 
built whose principal scheme is shown in Fig. HJ The basic element of this experimental set- 
up is computer data acquisition (DAQ) board NIPCI 6251 from National Instruments. This 
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board has two analog outputs and 16 analog inputs. The operating voltage is ±10 V, the 
maximal registration speed is 1 mus. The board has independent output and input buffer for 
4000 point. In our experiments, the output signal 10 V was not sufficient and the Wideband 
Power Amplifier KH model 7600 from Krohn-Hite Corporation with amplification coefficients 
5 and 25 times was used. It gives a possibility to have the output signal ±250 V. For the 
input signal this board has internal amplifier with coefficients 1, 2, 4, 8, 16, 32, 64, and 128 
times. A photo-detector was connected to input board plate for optical measurements. 

The software for experimental set-up has built-in functions for analog output-input. The 
program has three functional blocks that make operations with the set-up highly effective. 
The first block is a programmable generator that realizes any form of signal with duration 
of 2000 points. The duration of one point can be set from 1 fis to 1 s. The second block is 
a measuring block, which save 4000 values of the input voltage with step from 1 /is to 1 s. 
The operation of the first and the second blocks is synchronized inside the DAQ board and 
cannot be disturbed by computer interruptions. The third block is used to accumulate the 
experimental data during the working period. 



B. Theory 

In this section we apply our general theoretical approach described in Sec. [IT] to the 
special case of DHFLC polarization gratings. As is schematically shown in Figs. [1] and [31 in 
deformed helix ferroelectric (DHF) liquid crystals, the director 

d = (d x , d y , d z ) = cos 9 t x + sin 6 t cos $ y ± sin 6 t sin $ z (39) 

lies on the smectic cone with the helix tilt angle 6t and rotates in a helical fashion about a 
uniform in-plane twist axis (the x axis) forming the DHF helix. 

In the presence of weak electric field, E = E z, which is well below its critical value at 
the helix unwinding transition, E <C E c , the azimuthal angle around the cone, $, can be 



written in the form 30, 32 



2tt 

$ = + A$(0) ss + a E sm(j), = — x, (40) 

where the electric field induced distortion A$ linearly depend on E through the electric field 
parameter cie proportional to the ratio of the applied and critical electric fields: E/E c . 

From Eq. (I40p . it is clear that in the regime of weak electric field, the helix pitch P defines 
both the grating period, A g = P, and the grating wave number, k g = 2n/P. For the LC 
dielectric tensor ([I]) with the DHF helical configuration fl39l) . the anisotropy parameters that 
enter the matrices M mn whose block structure is described in Eq. (1181) are given by 



e±Vzz = ^— jj = ~i ] v = u a sin 2 9 t , (41) 

1 + u a d 2 z 1 + v sm z $ 

Paz = fiza = z ; ; t - ^' (42) 



1 - 1 - v sin <P 



According to the computational procedure developed in the previous section, after insert- 
ing the above parameters into the matrices of the system (fT6"j) for the Floquet harmonics, 
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we need to find the evolution operator by solving the initial value problem (l2"Uj) . Then 
computing the transmission and the reflection matrices (I30p involve the following steps: 
(a) evaluation of the linking matrix ([55]) ; (b) calculation of the efficiency matrices (1351) by 
solving the system of linear equation fl34|) : (c) computing the matrices T n and R„ from the 
formulas fl371) and fl38|) . respectively. 

As a consequence of the anisotropy induced mode coupling, there are an infinite number of 
coupled matrix equations in the system f|T6|) . So, the evolution operator cannot be generally 
computed in the closed form without resorting to the methods of the perturbation theory 
or numerical analysis. 

In this paper, we shall restrict ourselves to the case where the helix pitch is smaller than 
the wavelength, P < A, so as to interpret the experimental data measured in DHFLC cells 
with the subwavelength pitch. In our experiments, the incident light travels in the air with 
and normally impinges on the cell. At n m = 1 and k p = 0, the condition 

q 9 = k g /k v&c = \/P > 1 (44) 

implies that q n > n m . In this case, from the above discussion after Eq. (1291) . there are no 
diffracted waves propagating in the air and nonzero orders of the transmitted and reflected 
beams are evanescent. In the zero-order approximation, where M mn = £ m0 <5 n oM o, the 
evanescent waves are neglected and the DHFLC cell appears to be effectively described as 
a uniformly anisotropic layer characterized by the matrix 

"-<*-(8fc8fc). (45) 

where (...) = (27r) _1 / ... d0. The elements of the 2x2 block matrices (M)^- are given 

J o 

by 

<MC = <M>2? = -</?„>#\ (46) 

(M)^ } = - <?2° lf\ (47) 

<MC = - Ar 1 ^, pp = ix qpj (48) 

where k vac q p = k p and = (tfg) is the averaged effective dielectric tensor. 



1. Effective dielectric tensor 

Our next step is to evaluate the averages that enter the matrix (1451) . We note that, from 
Eq. (j42]) and the result 

(sin$ [1 + wsin 2 $] _1 ) = (sin $ cos $ [1 + wsin 2 $] _1 ) = 0, (49) 

the averages ((3 az ) are zero. When the averaged matrix fj45|) at ((3 az ) = is compared 
with the matrix of a uniformly anisotropic layer (see, e.g., Eq. (25) in Ref. [26]), where 

(Paz) = taz/tzz, (i]zz) = 1 / '^zz an d e^g = e^g , it immediately follows that the effective 
anisotropic medium is biaxial. In this medium, the diagonal element 



4f ] = (Vzz)- 



(50) 
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gives the principal value of the effective dielectric tensor for the optic axis normal to the cell 
(parallel to the z axis). 

For our purposes, it will suffice to expand the averages up to second order terms in the 
field parameter 

([l + usin 2 ^- 1 ) « [l + t;]- 1/2 (l+^ 2 «|), (51) 
(cos 2 $ [1 + v sin 2 « 7„(1 + [1 + if /2 7„a|), (52) 
(cos$ [1 + wsin 2 ^]- 1 ) « -7„a B , 7„ = (\/ITw - l)/v, (53) 

so that the in-plane components of the dielectric tensor are given by 

e xP = lx y a E , Ixy = ~ cos 9 t sin t 7„(e|| - e±), (54a) 
e£f > = c + Ae = e (°2 + 7 xx«|, (54b) 
e y f = e-Ae=e§j+ lyy a 2 E , (54c) 



where 



The relation 



4 U J = ex(l + K -«)[! + ^]" 1/2 ), eg = e ± [l + ^ i/2 , (55) 
7xx = e±(u a - v)v[l + v]~ 1/2 7 2 , j yy = e±v[l + v] 1/2 j 2 , (56) 



tan(20 d ) = -f xy a E /Ae (57) 



defines the azimuthal angle, <f>^ giving orientation of the in-plane optical axes. 

Similar to the LC dielectric tensor (PQ), the tensor f)54p can be expressed in terms of the 
"director" d eff = (a£ ,d y , 0) = (cos0 d , sin0 d , 0) and its eigenvalues 



e± = n ±/{i = e± Aey'l + tan 2 (20 d ) (58) 

as follows 

C = ^ + (e--e + )4 efr) 4 eff) - (59) 



2. Light transmittance 



For the case of normal incidence, the transmission and reflection matrices can be easily 



obtained from the results of Refs. [36|, |37J in the limit of the wave vectors with vanishing 



tangential component (see also Eq (1B47j) derived in Appendix IB 3j) . We shall need to write 
down the resultant expression for the transmission matrix 



t± 



^xx ^xy 
^yx ^yy 



1 — p± exp(2in±h) 



exp(m±/i), p± 



Rt(20 d ) • <r 3 , 

i: 

n±/n + n m /fi 



(60) 
(61) 
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When the incident wave is linearly polarized along the x axis (the helix axis), the transmit- 
tance coefficient 

T xy = M 2 = |t+ ~ t ~ |2 sin 2 (20 d ), sin 2 (20 d ) = ^ , (62) 

where h = k va , c D is the thickness parameter, describes the intensity of the light passing 
through crossed polarizers. Note that, under certain conditions such as |p±| ~ 1, the trans- 
mittance f )62|) can be approximated by simpler formula 

T xy « sin 2 (5/2) sin 2 (20 d ), (63) 

where 5 = An e gh = (n + — n_)h is the difference in optical path of the ordinary and 
extraordinary waves known as the phase retardation. 




Applied electric field (V/um) 



Figure 5. TVansmittance of light passing through crossed polarizers, T xy , as a function of the 
applied elecric field for the DHFLC cell of thickness D = 130 fim filled with the FLC mixture 
FLC-576A. Parameters of the mixture are: n Q = 1.5 (n e = 1.72) is the ordinary (extraordinary) 
refractive index and 8t = 32 deg is the tilt angle. The experimental points are marked by squares. 
Solid line represents the theoretical curve computed for the electric field parameter oe = ^eE 
with 7# Ri 0.62 fim/V. 



C. Experimental results and modeling 

Now we turn back to the electro-optical measurements in the DHFLC cells described in 
Sec. IIII Al The transmittance versus electric field curve shown in Fig. [5] presents the results 
measured at the wavelength of light generated by the He-Ne laser with A = 650 nm in the 
cell in which the thickness of the DHFLC layer, D, was about 130 /j,m. 
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The formula for the transmittance (16"2"j) can be combined with the expressions for the 
principal values of the effective refractive indices (j58p to evaluate dependence of T xy on the 
applied electric field, E. The known parameters characterizing the FLC mixture FLC-576A 
that enter our formulas are the ordinary (extraordinary) refractive index and the tilt angle 
estimated at n Q = 1.5 (n e = 1.72) and 6 t = 32 deg, respectively. So, for the anisotropy 
parameters, u a and v, we have: u a ~ 0.315 and v = u a sin 2 9 t w 0.09. 

From the discussion after Eq. fHOjl . the electric field parameter is proportional to 
the electric field, a# = JeE, and thus is determined by the coefficient of proportionality 
7e- This coefficient is the only fitting parameter in our calculations. Figure shows that 
the theoretical curve computed at 7^ ~ 0.62 /zm/V and the experimental data are in close 
agreement. 



0-14i | | | | | | | 




Applied voltage (V) 

Figure 6. Light transmittance, T xy , versus the voltage, U, applied across the DHFLC cell filled 
with the FLC mixture FLC-576A. The cell thickness is estimated at about D = 50 ± 2 [im. The 
experimental points are marked by squares. The theoretical curve (solid line) is computed at the 
field parameter cxe = "fvU with 7y ~ 0.0124 V . 

Similarly, in Fig. El the transmittance T xy is plotted against the applied voltage, U, for 
the case where the cell thickness is about 50 ± 2 /im. The theoretical results that give a 
good fit to the data of measurements are evaluated at the field parameter = jyU with 
7v ~ 0.0124 V -1 . So, in agreement with the above estimate, the cell thickness can be 
assessed at about D = "iE/lv ~ 50 /im. 

Referring to Figs and El the transmittance oscillates with the magnitude of the applied 
field. From Eqs. fIBTl) and fl62l) . this is the dependence of the effective refractive indices, n + 
and n_, on the electric field parameter, a^, that manifests itself in these oscillations. More 
precisely, as is evident from the approximate expression for the transmittance (163]) . they are 
due to variations in the phase retardation, 5, arising from the electric field induced birefrin- 
gence, An e g = n + — 72_. So, loci of extrema (minima and maxima) of the transmittance, 
T xy , can be used to obtain dependence of the birefringence on the applied voltage (more 
details on this method can be found, e.g., in Ref. (38j). 
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n 1 1 1 ,0.02 




Applied voltage (V) 

Figure 7. Electrically controlled birefringence, An e g(U) — An e ff(0), as a function of applied voltage. 
The theoretical curve (solid line) is computed from Eq. ([58]) at = j\/U and 7y 0.0124 V . 



Figure [7] shows the results for the voltage dependence of the electrically controlled bire- 
fringence, An e g(U) — An e ff(0), measured for the mixture FLC-576A. As can be seen from 
the figure, there is a remarkable accord between the experimental points the curve evaluated 
from the principal values of the effective dielectric tensor (15 9p . 

In addition to the birefringence, the transmittance ( 16 2 j) depends on the azimuthal angle, 
0d, given in Eq. ( ]5"Tj) . This angle describes the electric field induced rotation of the optical 
axes of the dielectric tensor (151?]) and, under the action of the electric field, E, its magnitude 
varies in the range between zero and 7r/4. 

A quick inspection of the formula (lrJ2"]) shows that the transmittance will always vanish 
in the zero-field limit with E = 0. Such behavior, however, comes as no surprise if we recall 
that the electric field of the incident wave is assumed to be parallel to the helix axis, Ei nc || x, 
which is the zero-voltage optical axis. 

Interestingly, the transmittance versus electric field curve will flatten in the vicinity of 
the origin when the cell thickness is chosen in such way that, in addition to sin0d, the 
birefringence dependent factor equals zero at E — 0. The latter is the case for the curves 
shown in Fig. [6] 



IV. DISCUSSION AND CONCLUSIONS 



In this paper we have formulated the theoretical approach to the optical properties of 
polarization gratings that can be regarded as an extension of the method developed in 
Refs IT 



26 



for stratified anisotropic media. Mathematically, in this approach, the key 
formulas (137]) and (l38j) give the transmission and reflection matrices of diffraction orders, 
T n and R n , expressed in terms of the linking matrix (155]) which is related to the evolution 
operator (|2"01 of the system of matrix equations ffTB"]) for the Floquet harmonics (|T4l de- 
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rived from the Maxwell's equations for the lateral components of the electric and magnetic 
fields (1T0]). 

This method goes beyond the well known limitations of the Jones matrix formalism at the 
expense of simplicity. To some extent it can be regarded as a version of the coupled mode 
analysis and our derivation of equations for the Floquet harmonics bear some similarity to the 



differential theory of light diffraction [39| which is based on the integration of a differential 
set of equations derived from Maxwell's equations projected onto some functional bases. 

In contrast to a theoretical method recently suggested in [Io[ for a reverse twisted nematic 
LC grating, our approach is formulated without recourse to the vector theory of scattering 
for a far field. Note that the optical properties of PGs were also analyzed numerically using 
the finite-difference time-domain method in Ref. 4l[. Although such analysis is general and 
useful, it can be very computationally intense so as to produce the results of sufficiently high 
accuracy. 

We have used the method described in Sec. [Ill as a tool of theoretical investigation into the 
electro-optic properties of deformed helix ferroelectric liquid crystal gratings with subwave- 
length pitch. It was previously demonstrated that a short helix pitch FLC (P = 0.4 -0.8 /im) 
provides an effective phase shift chan ge o f a transmitted polarized light beam as a function 
of the applied electric field intensity [27J. The new mixtures with the helix pitch being in 
UV region j35[ allow to get rid of complications related to the light scattering effects. 

So, in such mixtures, we deal with a pure electrically controlled phase shift plate based on 
the electro-optical mode called as deformed helix ferroelectric (DHF). DHF is a convenient 
operation mode able to ensure both low voltage and fast switching liquid crystalline light 
shutters (the response time is less than 200 /is at driving electric field around 1 V//xm). For 
these reasons, it is important to understand behavior of the electrically controlled birefrin- 
gence in such DHFLC cells. 

We have shown that, in the short-pitch approximation, DHFLC cells are equivalent to 
uniformly anisotropic biaxial layers with the optical axis normal to the substrate plane. 
The latter implies that normally incident light will feel only uniaxial in-plane anisotropy in 
agreement with the recent results on the ellipticity of light transmitted through a DHFLC 
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We have computed the averaged dielectric tensor as a function of the applied electric field 
and have used the results to evaluate the light transmittance measured in our experiments. 
A comparison between theoretical and experimental results presented in Figs. [5H7] shows 
that the predictions of the theory are in good agreement with the experimental data. Note 
that our calculations of the averaged dielectric tensor rely on the approximate expression for 
the azimuthal angle ( 140|) which is applicable for sufficiently low voltages. A more accurate 
description of the DHFLC orientational structure is required when the voltage increases 
approaching the unwinding transition [431 ]. 
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Appendix A: Derivation of equations for lateral components 

In this section we discuss how to exclude the ^-components of the electromagnetic field, 
E z and H z , that enter the representation (jHJ), from the Maxwell equations (j3J). Our task 
is to derive the closed system of equations for the lateral (tangential) components, Ep and 
Hp. 

We begin with substituting Eqs. (jSJ) and (jHJ) and have Maxwell's equations @ recast into 
the following form: 

- id T [z x E P ] = /iH — V p x E, (Ala) 

— id T Hp = D + V p x H, (Alb) 

where the explicit expressions for the last terms on the right hand side of the system flAlj) 
are as follows 

V p x E = -V^E Z + (V^ • Ep) z, (A2a) 
V p x H = -V^H Z + (V p • Hp) z. (A2b) 

We can now substitute the electric displacement field written as a sum of the normal and 
in-plane components 

D = D z z + Dp, (A3) 
into Eq. (I Alb)) and derive the following expression for its z-component 

D z = e zz E z + (e 2 ■ E P ) = - ( V p ■ H P ) , (A4) 

where e z = (e zx ,e zy ). 

From Eqs. flAlj) and flA4j) . it is not difficult to deduce the relations 

H z = f i- 1 (V^-E P ) (A5) 
^ = -e- 1 [(e z -Ep) + (Vp-Hp)] (A6) 

linking the normal (along the z axis) and the lateral (perpendicular to the z axis) compo- 
nents. 

By using the relation (jA6D . we obtain the tangential component of the field flA30 



Dp = e' z E z + e z -Ep = e P -E P -e' z e" 1 (V p • H P ) (A7) 

where e z — I €xx €xy J ; e' z — (e xz , e yz ) and the effective dielectric tensor, ep, for the lateral 

\ e yx e yy J 

components is given by 

e P = e z - e- z l e' z ®e z . (A8) 

Maxwell's equations AAID can now be combined with the relations (IA2I) to yield the 

system 

- id T E P = //Hp + Vp£ 2 , (A9a) 

-z9 r Hp = D P - X?iH z , (A9b) 
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where H z , E z and Dp are given in Eq. (1A5[) . Eq. (1A6j) and Eq. (IA7[) . respectively. 
So, this system immediately gives the final result 

- id T E P = -V p [e- Z l (e z ■ E P )] + yR P - Vpfe" 1 (V p - Hp)], (AlOa) 

- id T Hp = ep • Ep - V p ± [(V p ± ■ E P )//*] - e' z e" 1 (V p ■ H P ) (AlOb) 

that can be easily rewritten in the matrix form ( jiop used in Sec. [TT1 



Appendix B: Stratified anisotropic media 

In this appendix we the limiting case of stratified medium with vanishing modulation of 
the dielectric tensor (IT51) . so that e = Eq. Then the only mode of interest is the zero-order 
harmonics F = Fo governed by the matrix equation 

-i9 T F(r) = M(r).F(r), (Bl) 

where M = M 00 . 

The evolution operator of Eq. (IBip is defined as follows 

-id T U(r, r ) = M(r) • U(r, r ), U(r , r ) = I 4 . (B2) 

We now briefly touch on some of its basic properties. 



1. Operator of evolution 

We begin with the relation 

U(r,ro) = U(r,ri).U(ri > r ) (B3) 

known as the composition law. This result derives from the fact that the operator U(r, r ) • 
U _1 (ti, r ) is the solution of the initial value problem HB2[) with r replaced by 

From the composition law (1B3|) it immeadiately follow that the inverse of the evolution 
operator is given by 

U- 1 (r,r ) = U(r ,r) (B4) 

and can be found by solving the initial value problem 

idrU-^To) = U-^r.ro) • M(r), U-^ro.ro) = I 4 - (B5) 

For non-absorbing media with symmetric dielectric tensor, the matrix M is real-valued, 
M* = M, and meets the following symmetry identities 2o| : 



M* = M, (G ■ M) = G ■ M, G = I ^ , (B6) 

where an asterisk and the superscript T indicate complex conjugation and matrix trans- 
position, respectively. In this case, the evolution operator and its inverse are related as 
follows: 

U- 1 (r,r ) = G-U t (r,r )-G, (B7) 
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where a dagger will denote Hermitian conjugation. By using the relations (1B6[) . it is not 
difficult to verify that the operator on the right hand side of Eq. (IB7I) is the solution of the 
Cauchy problem (1B5p . 

We are now in position to deduce the unitarity relation 



W 1 = G 3 W f G 3 



wJi -w 



wU w 



21 

t I ' 



12 T, 22 



G 3 = diag(I 2 , -I 2 ) 



(B8) 



for the linking matrix 



w = v- 1 -u- 1 (/i,o).v 



(B9) 



by using Eq. flBTI) in combination with the algebraic identity 

V™ • G • V m = iV m G 3 , 



(BIO) 



where iV m = 2g m /yU m , for the eigenvector matrix 



V m = V m (g n 



E m — cr 3 E m 
1 1 cr 3 H m 



(Bll) 



with E m = diag(g m /n m , 1) and /i m H m = diag( 



z=0 




£ ij s J-i ij 



z=D 



^ m> Mr 




- Z 




Figure 8. Four- wave geometry with two incident waves, E-j^ and E^J, impinging onto the entrance 
(z = 0) and exit [z = D) faces, respectively. 
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2. Linking Matrix 

For closer examination of properties of the linking matrix (IB9j) . we consider the four- wave 
geometry shown in Fig. El where k^? = k tr = k inc = k trm and k^J = k rt = k re fl. In this 
geometry, there are two plane waves, and E[~ c \ incident on both bounding surfaces of 
the anisotropic layer. For the transmitted and reflected waves, E tr and E rt , excited by the 
incident light, the input-output relation (l30i) takes the following generalized form: 



E tr 

E rt J \ R i 




[B12) 



where T + (R+) is the transmission (reflection) matrix defined in Eq. (1301) . whereas the 
mirror symmetric case where the incident wave is impinging onto the exit face of the sample 
is described by the transmission (reflection) matrix T_ (R ). 
The matrix (1B9I) is introduced through the relation 25, 26| 



W ( Jfr) = m (B13) 



E(- J \ T? / 

inc / V ^rt / 

linking the electric field vector amplitudes of the waves in the half spaces z < and z > D 
bounded by the faces of the layer. From Eqs. (1B12I) and (IB13I) . its block structure expressed 
in terms of the transmission and reflection matrices is as follows 

w =(r + T Vt_-r;'t^.r_)' < B14 > 

Similarly, for inverse of the linking matrix, we have 

, _ / w«7" w<^>\ _ (t + - R_ ■ t: 1 ■ r + r_ ■ t:'\ 
" w<7"w«-> -{ -t:'.r + t:' )~ < B15 > 



We can now use the unitarity relation (1B8j) for non-absorbing layers to derive the energy 
conservation law 

Tl • T± + R4 • R± = I 2 . (B16) 
along with the relations for the block matrices 

W n = T;\ W 22 = [TlY, (B17a) 
Wia = -T; 1 ■ R_ = [TI 1 ■ R + ]\ (B17b) 

w 21 = r + ■ t; 1 = -[r_ ■ t: 1 ] 1 . (Bi7c) 

In the translation invariant case of uniform anisotropy, the matrix M is independent of 
r and the operator of evolution is given by 

U(r, r Q ) = U(r - r ) = exp{*M (r - r )}. (B18) 
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Then, the unitarity condition 26 



IT 1 = U*, W" 1 = W* (B19) 

can be combined with Eq. ( \B8\i to yield the additional symmetry relations for Wy 

W? = W«, Wf 2 = -W 21 (B20) 

that give the following algebraic identities for the transmission and reflection matrices: 

= T±, R+ = R_, (B21) 
= -R* T ■ T T • R" 1 . (B22) 

3. Uniformly anisotropic planar structure 

In this section we present the results for the planar oriented nematic cell. For the planar 
orientational structure, the director is given by 

d = (d x , d y , d z ) = (cos d , sin <f> d , 0), (B23) 

so that the layer is uniformly anisotropic. 

Assuming that the plane of incidence is defined as the x-z plane (0i nc = 0), we write 
down the elements of the matrix M as follows 

M 12 =^I 2 -^4l 2 + o- 3 ), M ii = 0, (B24) 

2 

M 2 i = -f^s - o-a) + e c {l 2 + u a [cos(20 d ) <r 3 + sin(20 d ) ct x ] }, (B25) 
e c = (e„+e ± )/2, u a = q x = q x p \ (B26) 

where the azimuthal angle of the in-plane optical axis </> d also gives the angle between the 
director and the plane of incidence. 

For uniform anisotropy, this matrix is constant. Therefore, the operator of evolution can 
be expressed in terms of the eigenvalue and eigenvector matrices, A and V, as follows 

U(h) = exp{zM h} = \ ■ exp{i\ h} ■ V" 1 . (B27) 

It is not difficult to solve the eigenvalue problem for the matrix M and find the expressions 
for the eigenvalues that enter the eigenvalue matrix 

A = diag(Q, -Q), Q = diag(g e , q ), (B28) 



y/n*-ql(l + u a dl), q = Vn 2 o-Ql (B29) 

where u a = Ae/e±. Similarly, after computing the eigenvectors, we obtain the eigenvector 
matrix in the following form: 



E = ^ Ml - - gl/nl) d q \ (B31) 

dxQe dyTl 

d y q e -d x [n 2 -q 2 x ] 



H =fe e A d S l , 21 ). (D32) 
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Upon substituting Eqs. (1B27[) — f]B32|) into Eq. (jB9j), some rather straightforward algebraic 
manipulations give the linking matrix 



W = iV m 1 diag(I 2 , <r 3 ) ■ W • diag(I 2 , ct z 



W 

w d 



A T + -A T 
-A T AT 



A+ A 

A- A +/ 

W 

o w +/ 

A± = E m • H ± H m • E, N = diag(iV e , N ) 



(^o~ w + ) ' w± = e M±*Qh] ■ n- 1 , 



q x d x ), N = 2q fi(n - q x d x ) 



(B33) 
(B34) 

(B35) 
(B36) 
(B37) 



where N m = 2q m / fi m . From Eq. (1B34j) . the block 2x2 matrices of W are given by 

Wh = N m W n = A + ■ W • A^ - A_ • W + ■ AZ, 

W 22 = N m a 3 ■ W 22 ■ 0-3 = A + ■ W + ■ A^ - A_ • W_ ■ A T , 

W 21 = N m cr 3 ■ W 21 = -W T 2 = -NjWu ■ <rz] T = 



A_ • W_ • AT - A + • W+ • A J 



(B38a) 
(B38b) 

(B38c) 



Finally, we can combine Eq. (1B14I) with Eq. (1B33|) to derive the expressions for the 
transmission and reflection matrices 



A^wr/, r+ = <r 3 W21 ■ wr/. 



(B39) 



As it can be seen from the formulas (IB38[) . the symmetry relations (IB2Q[) are satisfied. 
Interestingly, when the eigenvalue matrix (IB28I) is real so that = W_ and A^_ = A^, 
close inspection of the expressions (IB38I) shows that, in addition to the symmetry relations 
for uniform anisotropy flB20j) . the case of uniform planar structure is characterized by the 
following algebraic identities: 

W| 2 = <r 3 • W n ■ 0-3, Wj x = -0-3 • W 21 • <t 3 . (B40) 

These identities and the unitarity conditions (IB 171) can now be used to deduce the relations 

T+ = 0-3 • T_ • tr 3 , R+ = cr 3 R_ cr 3 (B41) 

linking the transmission (reflection) matrix, T + = T (R + = R), and its mirror symmetric 
counterpart T_ (R_). 

Before closing this section we briefly comment on the important special case of normal 
incidence that occurs at q x = 0. In this case, the matrices A± defined in Eq. flB36j) can be 
written in the factorized form 



A ± (0 d )=Rt(0 d )-A ± (O)= ^ fA ■ 



—n r 





/i m n Q ± n n n 



\ 



(B42) 
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cos cb sin 

where Rtf(f>) = I . , 7 1 is the matrix describing rotation about the z axis 

1 sin0 cosff ' 



angle 0. Substituting Eq. (1B42|) into Eq. (1B38|) gives the block matrices 



W 



ijWd 



) = Rt( 



Wy (0) • Rt(-0 d ) 



expressed as a function of the director azimuthal angle d . 
The result for the transmission and reflection matrices 

T ± (0 d ) = Rt(±0 d ) • T(0) • Rt(=F0 d ), R±(0d) = Rt(=F0d) • R(0) • Rt(=F&) 

where the diagonal matrices T(0) = diag(t e ,t Q ) and R(0) = diag(r e , — r Q ) 

1 — Pa ,. ,n 1 — exp(2in a /i) 



1 - p 2 a exp(2in a h) 
n a /n - n m /p 
n a /n + n m /p 

in 



exp(in a h), r 
a E {e, o] 



1 - p 2 a exp(2in a h) 



Pa, 



describe the case in which the director flB23j) lies in the incidence plane, immediately 
from the relations f lB39j) and flB41j) . Finally, the expressions for the matrices (IB44I) 



T±(0 d ) 

R±(0 d ) 



te + ^0 
2 

r e + r 



h + Rt(±20 d ) • a 3 , 



o"3 + ! ^^Rt( T 20 d ) 



2 J 2 

can be readily obtained by using the identity: cr 3 ■ Rt(0) ■ cr 3 = Rt(- 



by the 
(B43) 

(B44) 

(B45) 
(B46) 
follows 

(B47) 
(B48) 
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